Anatomy of the western Java plate interface from depth-migrated seismic images

Newly pre-stack depth-migrated seismic images resolve the structural details of the western Java forearc and plate interface. The structural segmentation of the forearc into discrete mechanical domains correlates with distinct deformation styles. Approximately 2/3 of the trench sediment fill is detached and incorporated into frontal prism imbricates, while the floor sequence is underthrust beneath the décollement. Western Java, however, differs markedly from margins such as Nankai or Barbados, where a uniform, continuous décollement reflector has been imaged. In our study area, the plate interface reveals a spatially irregular, nonlinear pattern characterized by the morphological relief of subducted seamounts and thicker than average patches of underthrust sediment. The underthrust sediment is associated with a low velocity zone as determined from wide-angle data. Active underplating is not resolved, but likely contributes to the uplift of the large bivergent wedge that constitutes the forearc high. Our profile is located 100 km west of the 2006 Java tsunami earthquake. The heterogeneous décollement zone regulates the friction behavior of the shallow subduction environment where the earthquake occurred. The alternating pattern of enhanced frictional contact zones associated with oceanic basement relief and weak material patches of underthrust sediment influences seismic coupling and possibly contributed to the heterogeneous slip distribution. Our seismic images resolve a steeply dipping splay fault, which originates at the décollement and terminates at the sea floor and which potentially contributes to tsunami generation during co-seismic activity

Movement along a low-angle normal fault: The S reflector west of Spain

[1] The existence of normal faults that moved at low angles (less than 20°) has long been debated. One possible low-angle fault is the S detachment at the west Galicia (Spain) margin and thought to occur at the top of serpentinized mantle. It is unlikely that S was a large submarine slide as it was probably active over several million years without the development of any compressional features such as toe thrusts, it appears to have rooted beneath the conjugate Flemish Cap margin, and it is similar to structures elsewhere that also appear to be rooted detachments. Here we analyze depth images to identify synrift sediment packages above S and use the geometry of these synrift packages to constrain the angle at which S both formed and remained active. We find that S must have remained active at angles below 15°, too low to be explained simply by the low friction coefficient of partially serpentinized peridotites. Instead, we suggest that transient high fluid pressures must have developed within the serpentinites and propose a model in which anastomosing fault strands are alternately active and sealed, enabling moderately high fluid pressures to develop

Marine forearc structure of eastern Java and its role in the 1994 Java tsunami earthquake

We resolve a previously unrecognized shallow subducting seamount from a re-processed multichannel seismic depth image crossing the 1994 M7.8 Java tsunami earthquake slip area. Seamount subduction is related to the uplift of the overriding plate by lateral shortening and vertical thickening, causing pronounced back-thrusting at the landward slope of the forearc high and the formation of splay faults branching off the landward flank of the subducting seamount. The location of the seamount in relation to the 1994 earthquake hypocentre and its co-seismic slip model suggests that the seamount acted as a seismic barrier to the up-dip co-seismic rupture propagation of this moderate size earthquake. The wrapping of the co-seismic slip contours around the seamount indicates that it diverted rupture propagation, documenting the control of forearc structures on seismic rupture

Stochastic Heterogeneity Mapping around a Mediterranean salt lens

We present the first application of Stochastic Heterogeneity Mapping based on the band-limited von Kármán function to a seismic reflection stack of a Mediterranean water eddy (meddy), a large salt lens of Mediterranean water. This process extracts two stochastic parameters directly from the reflectivity field of the seismic data: the Hurst number, which ranges from 0 to 1, and the correlation length (scale length). Lower Hurst numbers represent a richer range of high wavenumbers and correspond to a broader range of heterogeneity in reflection events. The Hurst number estimate for the top of the meddy (0.39) compares well with recent theoretical work, which required values between 0.25 and 0.5 to model internal wave surfaces in open ocean conditions based on simulating a Garrett-Munk spectrum (GM76) slope of −2. The scale lengths obtained do not fit as well to seismic reflection events as those used in other studies to model internal waves. We suggest two explanations for this discrepancy: (1) due to the fact that the stochastic parameters are derived from the reflectivity field rather than the impedance field the estimated scale lengths may be underestimated, as has been reported; and (2) because the meddy seismic image is a two-dimensional slice of a complex and dynamic three-dimensional object, the derived scale lengths are biased to the direction of flow. Nonetheless, varying stochastic parameters, which correspond to different spectral slopes in the Garrett-Munk spectrum (horizontal wavenumber spectrum), can provide an estimate of different internal wave scales from seismic data alone. We hence introduce Stochastic Heterogeneity Mapping as a novel tool in physical oceanography

Resolving the fine-scale velocity structure of continental hyperextension at the Deep Galicia Margin using full-waveform inversion

Continental hyperextension during magma-poor rifting at the Deep Galicia Margin is characterised by a complex pattern of faulting, thin continental fault blocks, and the serpentinisation, with local exhumation, of mantle peridotites along the S-reflector, interpreted as a detachment surface. In order to understand fully the evolution of these features, it is important to image seismically the structure and to model the velocity structure to the greatest resolution possible. Travel-time tomography models have revealed the long-wavelength velocity structure of this hyperextended domain, but are often insufficient to match accurately the short-wavelength structure observed in reflection seismic imaging. Here we demonstrate the application of two-dimensional (2D) time-domain acoustic full-waveform inversion to deep water seismic data collected at the Deep Galicia Margin, in order to attain a high resolution velocity model of continental hyperextension. We have used several quality assurance procedures to assess the velocity model, including comparison of the observed and modelled waveforms, checkerboard tests, testing of parameter and inversion strategy, and comparison with the migrated reflection image. Our final model exhibits an increase in the resolution of subsurface velocities, with particular improvement observed in the westernmost continental fault blocks, with a clear rotation of the velocity field to match steeply dipping reflectors. Across the S-reflector there is a sharpening in the velocity contrast, with lower velocities beneath S indicative of preferential mantle serpentinisation. This study supports the hypothesis that normal faulting acts to hydrate the upper mantle peridotite, observed as a systematic decrease in seismic velocities, consistent with increased serpentinisation. Our results confirm the feasibility of applying the full-waveform inversion method to sparse, deep water crustal datasets

Deep structure of the Ionian Sea and Sicily Dionysus - Cruise No. M111, October 10 - November 1, 2014, Catania (Italy) – Catania (Italy)

Summary The origin of the Ionian Sea lithosphere and the deep structure of its margins remain a little investigated part of the Mediterranean Sea. To shed light on the plate tectonic setting in this central part of southern Europe, R/V METEOR cruise M111 set out to acquire deep penetrating seismic data in the Ionian Sea. M111 formed the core of an amphibious investigation covering the Ionian Sea and island of Sicily. A total of 153 OBS/OBH deployments using French and German instruments were successfully carried out, in addition to 12 land stations installed on Sicily, which recorded the offshore air gun shots. The aim of this onshore-offshore study is to quantify the deep geometry and architecture of the Calabria subduction zone and Ionian Sea lithosphere and to shed light on the nature of the Ionian Sea crust (oceanic crust vs. thinned continental crust). Investigating the structure of the Ionian crust and lithospheric mantle will contribute to unravel the unknown ocean-continent transition and Tethys margin. Analyzing the tectonic activity and active deformation zones is essential for understanding the subduction processes that underlie the neotectonics of the Calabrian subduction zone and earthquake hazard of the Calabria/Sicily region, especially in the vicinity of local decoupling zones

Ionian Abyssal Plain: a window into the Tethys oceanic lithosphere

The nature of the Ionian Sea crust has been the subject of scientific debate for more than 30 years, mainly because seismic imaging of the deep crust and upper mantle of the Ionian Abyssal Plain (IAP) has not been conclusive to date. The IAP is sandwiched between the Calabrian and Hellenic subduction zones in the central Mediterranean. A NNE–SSW-oriented 131&thinsp;km long seismic refraction and wide-angle reflection profile, consisting of eight ocean bottom seismometers and hydrophones, was acquired in 2014. The profile was designed to univocally confirm the proposed oceanic nature of the IAP crust as a remnant of the Tethys and to confute its interpretation as a strongly thinned part of the African continental crust. A P-wave velocity model developed from travel-time forward modelling is refined by gravimetric data and synthetic modelling of the seismic data. A roughly 6–7&thinsp;km thick crust with velocities ranging from 5.1 to 7.2&thinsp;km&thinsp;s−1, top to bottom, can be traced throughout the IAP. In the vicinity of the Medina seamounts at the southern IAP boundary, the crust thickens to about 9&thinsp;km and seismic velocities decrease to 6.8&thinsp;km&thinsp;s−1 at the crust–mantle boundary. The seismic velocity distribution and depth of the crust–mantle boundary in the IAP document its oceanic nature and support the interpretation of the IAP as a remnant of the Tethys lithosphere with the Malta Escarpment as a transform margin and a Tethys opening in the NNW–SSE direction.</p

Resolving the fine-scale velocity structure of continental hyperextension at the Deep Galicia Margin using full-waveform inversion

Continental hyperextension during magma-poor rifting at the Deep Galicia Margin is characterized by a complex pattern of faulting, thin continental fault blocks and the serpentinization, with local exhumation, of mantle peridotites along the S-reflector, interpreted as a detachment surface. In order to understand fully the evolution of these features, it is important to image seismically the structure and to model the velocity structure to the greatest resolution possible. Traveltime tomography models have revealed the long-wavelength velocity structure of this hyperextended domain, but are often insufficient to match accurately the short-wavelength structure observed in reflection seismic imaging. Here, we demonstrate the application of 2-D time-domain acoustic full-waveform inversion (FWI) to deep-water seismic data collected at the Deep Galicia Margin, in order to attain a high-resolution velocity model of continental hyperextension. We have used several quality assurance procedures to assess the velocity model, including comparison of the observed and modeled waveforms, checkerboard tests, testing of parameter and inversion strategy and comparison with the migrated reflection image. Our final model exhibits an increase in the resolution of subsurface velocities, with particular improvement observed in the westernmost continental fault blocks, with a clear rotation of the velocity field to match steeply dipping reflectors. Across the S-reflector, there is a sharpening in the velocity contrast, with lower velocities beneath S indicative of preferential mantle serpentinization. This study supports the hypothesis that normal faulting acts to hydrate the upper-mantle peridotite, observed as a systematic decrease in seismic velocities, consistent with increased serpentinization. Our results confirm the feasibility of applying the FWI method to sparse, deep-water crustal data sets

Continental hyperextension, mantle exhumation, and thin oceanic crust at the continent-ocean transition, West Iberia: New insights from wide-angle seismic

Hyperextension of continental crust at the Deep Galicia rifted margin in the North Atlantic has been accommodated by the rotation of continental fault blocks, which are underlain by the S reflector, an interpreted detachment fault, along which exhumed and serpentinized mantle peridotite is observed. West of these features, the enigmatic Peridotite Ridge has been inferred to delimit the western extent of the continent‐ocean transition. An outstanding question at this margin is where oceanic crust begins, with little existing data to constrain this boundary and a lack of clear seafloor spreading magnetic anomalies. Here we present results from a 160 km long wide‐angle seismic profile (Western Extension 1). Travel time tomography models of the crustal compressional velocity structure reveal highly thinned and rotated crustal blocks separated from the underlying mantle by the S reflector. The S reflector correlates with the 6.0–7.0 km s−1 velocity contours, corresponding to peridotite serpentinization of 60–30%, respectively. West of the Peridotite Ridge, shallow and sparse Moho reflections indicate the earliest formation of an anomalously thin oceanic crustal layer, which increases in thickness from ~0.5 km at ~20 km west of the Peridotite Ridge to ~1.5 km, 35 km further west. P wave velocities increase smoothly and rapidly below top basement, to a depth of 2.8–3.5 km, with an average velocity gradient of 1.0 s−1. Below this, velocities slowly increase toward typical mantle velocities. Such a downward increase into mantle velocities is interpreted as decreasing serpentinization of mantle rock with depth